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Abstract 

In the high spin limit the minimal anomalous dimension of (fixed) twist operators in the si (2) 
sector of planar N = 4 Super Yang-Mills theory expands as j(g,s,L) = /(g) Ins + f s i(g,L) + 

oo 

S 7^(fiS-^) (his)"™ + • • •• We find that the sub-logarithmic contribution ^ n \g,L) is governed 

n=l 

by a linear integral equation, depending on the solution of the linear integral equations appearing 
at the steps n' < n — 3. We work out this recursive procedure and determine explicitly "f^ n \g, L) (in 
particular j^(g,L) = and 7( n )(g,2) = y( n \g,3) = 0). Furthermore, we connect the ^ n \g,L) 
(for finite L) to the generalised scaling functions, fn{g), appearing in the limit of large twist 
L ~ Ins. Finally, we provide the first orders of weak and strong coupling for the first ^ n \g,L) 
(and hence fn\g)). 
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1 A quick outlook 



The calculation of anomalous dimensions in the planar M = 4 Super Yang-Mills (SYM) theory re- 
ceived a considerable boost as a consequence of the hints for an underlying integrable structure (cf. for 
instance [U El El HI E])- This meant the identification of the Bethe(-Yang) Ansatz equations valid in 
the asymptotic regime of very long operators. The energies (i.e. the integrable asymptotic spectrum) 
corresponds to planar anomalous dimensions of single trace operators, when the number of the com- 
pounds is very large. This approximation would correspond to the on-shell (or IR) description in the 
usual relativistic field theories [6]. Interestingly, in this regime it was possible to predict, in a relatively 
simple fashion, some exact weak-coupling expansions of anomalous dimensions up to a certain number 
of loops; these were successfully compared with results of field theoretical perturbation theory (e.g. 
[HUHE] and references therein). In some lucky cases, this approximation gave also the impressive access 
to some exact results in the strong coupling regime, dominated by the dual string theory (cf. e.g. [9] 
and references therein). More in general, an anomalous dimension is affected by corrections coming 
from the finite size of a composite operator, but this on-shell approach has also furnished a basis for 
the off-shell (any length at any coupling) prediction of the recent Thermodynamic Bethe Ansatz [10J. 
All these integrability based results are precision tests for the more general AdS/CFT correspondence 

EH). 

In this paper we focus on the study of the twist si (2) sector of Af = 4 planar SYM whose operators 
have the form 

Ti(V s Z L ) + . . . , (1.1) 

where the s symmetrised, traceless, covariant derivatives T> act in all the possible ways on the L bosonic 
fields Z. s is just the Lorentz spin and L an su(4) R-charge equal to the twist (=(classical dimension) - 
s) . Peculiarly, the minimal anomalous dimension of (11. II) shows a leading high spin (fixed L) behaviour 
with the logarithm of s, having as coefficient the so-called universal scaling function, /(<?)□ [T2" | [T51 [T4] . 
Impressively, it was obtained from the solution of a linear integral equation directly derived from the 
asymptotic Bethe Ansatz (ABA) via the root density approach j5j. Moreover, it was carefully studied 
and tested both in the weak [U E] and strong coupling limit [151 ESI EH EE] • 

In the following, we wish to prove that the high spin expansion goes on as a series of logarithmic 
(inverse) power 

oo 

j(g, s, L) = f(g) lns + f sl (g, L) + X> (r %, L) (Ins)-" + O ((his)" 00 ) , (1.2) 

n=l 

i.e. it is a genuine large size expansion in the size parameter Ins. In this context the sub-leading 
(constant) contribution f s i(g,L) received already much attention. In [Tj5] it was shown to come from 
the solution of a non-linear integral equation (NLIE). Then, in [20] it was obtained starting from a 
linear integral equation (LIE). Explicit weak and strong coupling expansions are present in [21] and 



1 It does not depend on L and equals twice the cusp anomalous dimension of light-like Wilson loop. 
2 With O ((In s)~°°) we indicate terms going to zero faster than any inverse powers of Ins. 



2 



agree with string theory computations [22]. Importantly, it is believed that both f(g) and f s i(g,L) are 
exactly given by this approach based on the ABA without finite size corrections. In fact, for instance the 
findings of [23] have proved that, at least at twist two and up to four loops, wrapping corrections start 
contributing at order O ((lns) 2 /s). More in general, we may expect that all the logarithmic sub-leading 
contributions in (11.21) . 7^(<7, L), are exact as they can be elaborated from the ABA. 

In this letter we focus on ^ n \g,L). We prove that each j^(g,L) is still determined by a linear 
integral equation (for the density of roots and holes), which can be equivalently rewritten as an infinite 
dimensional (matrix) linear system. Things can be arranged in such a way that for any n the systems 
have the same kernel, differing only in the inhomogeneous terms. These forcing terms are expressed 
as linear combinations of those driving the systems for the "reduced coefficients" , defined by equations 
(4.23,24) of [21] and appearing in the study of the so-called generalised scaling functions f m (g) [T9| [24] 
and fm (g), describing the anomalous dimension 

oo oo oo 

7 (<7, 3, L) = In s ]T f n (g)j n + £(ln s)~ r ]T /« (g)f + O ((In s )"») , (1.3) 

n=0 r=0 n=0 

in the limit [T5] 

L — 2 

s — > oo , L — > oo , j = — ■ fixed. (1.4) 

In s 

In this way a recursive structure is set up such that the driving term in the system n for 7^(gf, L) 
involves the solutions to those with nl < n — 3. Therefore, it is possible to push the computations at 
any desired order n and give a general expression for the form of 7* n * ) (g, L). From the latter the higher 
generalised scaling functions fm\g) can be extracted @. 



2 All- loops ABA and the (N)LIE at high spin 

In a series of papers we proposed the NLIE [25J technique for the study of asymptotic Bethe Ansatz 
equations describing the AdS/CFT correspondence. An interesting case where such a technique was 
recently applied [20] is the calculation of anomalous dimensions of single-trace operators belonging to 
the si (2) sector of the M = 4 SYM theory. In such a specific case, if one looks at the high spin limit at 
fixed twist, one has two crucial simplifications. Firstly, the non-linear integral contributions to the NLIE 
are depressed and contribute with terms O ((Ins) -00 ) [191 [20]. Secondly - as stated in the introduction 
- wrapping effects are also depressed and start contributing at order 0((lns) 2 /s) [23]. This means that 
in order to work out all the sub-leading terms behaving as (lns)~ n , with n > 0, we can safely rely on 
the linear integral terms of the NLIE coming from the asymptotic Bethe Ansatz. 

In this letter we focus on the minimal anomalous dimension state, with spin s and twist L. Such a 
state is characterised by s Bethe roots, localised in an interval [—6, b] of the real axis and by L real 'holes'. 
The positions of both the roots and the holes are symmetric with respect to the origin. For what concerns 
the holes, two lie outside [—b, b], the remaining L — 2 concentrate near the origin, with no roots lying in 

3 This is indeed the BES kernel determining f(g) [5] and f s i(g,L) [20]. [2Tj . 
4 They can also be computed independently as in Appendix B 
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between. The positions Vk of both roots and holes satisfy the condition Z(vk) = 7r(2fc+l), with k integer, 
where Z(u) = —Z(—u) is the counting function; the positions of the internal holes Uh are unknowns in 
the NLIE and are determined by the (non linear) relations Z(v,h) = n(2h + 1 — L) , h = 1, . . . , L — 2 , 
which supplement the (linear integral part of the) NLIE. 

We move from upshots of [20]. The crucial point is that the forcing term F(u) defined in [20] 
approximates the counting function Z(u) in the large s limit, if we neglect terms of order O ((Ins) -00 ). 
Hence, the density of roots and holes a{u) = 4-Z{u) can be computed (up to this order) by using 
the linear integral equation (3.52, 4.10 of [20]) for F(u). In specific, splitting the forcing term in one 
and higher than one loops contributions, F(u) = F (u) + F H (u), the linear integral equation for the 
higher- loop part F H (u) ia§ 



u 



—iL In 



1 + 



2x~(u) 2 



2x+{u) 2 

+ i9(u, + i arctan(w 

7T 



L-2 /i 
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i arctan(w 
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r+oo 
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/ dv 


/ — oo 





2x+(m)x (v) 
2x~ (u)x+ (v) 



■ & nr \ 

+ i—v(u,v) 
dv 



[F (v) + F H (v)}+O((\ns)-° o ) . 



This equation comes from (4.10) of [20], after removing all the O ((Ins) °°) terms. It has to be supple- 
mented with the expression for F (u) - (3.52) of [20] - which in Fourier transform reads: 



ikF (k) 



k — e 2 cos(/cs/a/2) ^— \ 
-4vr^ hA + 2tt 2^ 



e iku h 



(0) 



_1M 
e 2 



2sirihiM ^ 2sinh- 2 



4vrln2(5(fc) + ((Ins)" 00 ) . (2.2) 



In d2U E2} uf ] represent the one loop contribution to the (^-depending position of the internal holes 
Uh- Now, since the function F{u) approximates the counting function Z{u) up to O ((lns)~°°) terms, 
in this approximation the positions of the internal holes Uh are determined from the conditions 



F(u h ) = ir(2h + 1 - L) , h = 1, 



2. 



(2.3) 



The key point which allows to understand the origin of the logarithmic terms (lns) _n is related to the 
behaviour of the position of the holes as a function of the spin s in the s — > oo limit. In such a limit 
for fixed u F(u) diverges logarithmically (and expands in (inverse) powers of Ins) as 



F(u) = F {n) (u)(\ns)- n + O ((Ins)" 00 ) 



(2.4) 



hence it is natural to suppose that - in order to fulfil (12.31) - the position of the holes has to expand in 
inverse powers of Ins: 



Uh 



^a n ,,(lns)^ + 0((lns)-° o ) 



(2.5) 



n=i 



5 For notations we refer to seminal papers [31 HI E] - 
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A systematic way to perform such an expansion order by order in a recursive way is given by the 
Faa di Bruno formula for the derivatives of a composite function. We already used this trick in [23] 
and we refer to that paper for technical details. Introducing the derivatives in zero of the function 
a(u) = cr(—u) = 4-F(u) and developing them in powers of Ins, 

^a{u = 0) = J>W(ln S )-", (2.6) 

n=—l 

(ai n ^ = when r is odd) the condition (12.31) for the holes eventually gives 

00 p+l p-r+2 , « • 

n(2h+l-L) = ^WE^'E^ E II ^T 1 ^ 

P=l r=l 0"i.-Jp-r+a} m=1 3m ' 

OO P— 1 P—I P—I i+1 / 

+ swEEA £ II { -^fr z < ("> 

p=l i=0 r=l {ii,-,J P -r-!+l} ™ =1 • 7m ' 

where the j m contained in the second term of the r.h.s. are constrained by the conditions Y^m=i 2 3 m = r ' 
YfmLt 2 ™!™ = P + 1, the ones in the third term by Y^Z\ l+1 3m = r, J2m=i l+1 m 3m =p — l Equating 
l.h.s. and r.h.s. at all orders in Ins we obtain the following recursive equation 

V (-1) p—r+l , p-l p-l (l) p—r—l+l , \ a 

a P+l , h = -E-rt E 11 -^-LL— E 11 — p— • 

r=l °0 {ii,...j p _H-i} m=1 '=0 r=l °0 {ji.-^-r-1+i} m=1 

ai.ft = tzX) , (2.8) 

^0 

where now the j m contained in the first term of the r.h.s. are constrained by the conditions Yl^Ii 1 3 m = 
r + !> Z]m~=Ti +1 = P + 1 and the ones in the second term by Y!m=\ 3m = r, YXZi +1 m 3m =p-l 
Equation (12.81) is solved by iterations, allowing to express the coefficients a n ^ in terms of the derivatives 
in zero of the density of roots and holes. The first three of them are 

ix{2h + 1 - L) ix{2h + 1 - L)ctJ 0) 



CT o l°o 
tt(2/i + 1-L) / 4 1} / 4 0) V n\ , rN2 



3 Linear integral equations for the logarithmic terms 

After we clarified the behaviour of the position of the holes for high s, we continue with the standard 
treatment of equation (12.11) . in order to study the contributions to the anomalous dimension proportional 
to (lns) _n . For simplicity's sake, from now on we do not write the O ((Ins) -00 ) terms. 
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Since the anomalous dimension at high spin is given [7] by j(g, s,L) = lim CTjf ( fc ) ; where (Jh(u 



k^O 



4-F (u), we write (12 .ip in terms of the Fourier transform anik), restrict to k > and introduce 



S(k) 



sinh 
ir\k\ 



\k\ L-2 



■\&H{k) — ii rrr [cos kuh — cos kuf^ } =>- j(g, s, L) — 2 lim S(k) . (3.1) 

I sinh — ' L -I J 

■ ■ blllu 2 h=i 

The function ( 13. ip satisfies the linear integral equation 

t r+oo 1, 

S(k) = -[1 - Jo(^fc)] - <? 2 / -e^k(V2gk, y/2gt) ■ 
& Jo n 

cos 

(3.2) 



r 7r/^ 1 — p 2 



1 - e 2 cos ^= 



+ 7T— ^tJZ i costu h ~ 1] }= Vlns K{V2gk,0) + Ag 2 J - — -K*(V2gk, V2gt) + 



2 h=i 



+ -[1- J (v / 2pfc)]+ V7 B ^:(v / 2pfc,0) + g 2 (L-2) / dte-^K(V2gk,V2gt)- 



1 — 



sinh I 



<7 2 /" + °° dtK(V2gk, V2gt) ^ Lh = l[c ° S ^ h 1] - <? 2 /" + °° dte^K(V2gk, V2 
Jo smh § ,/ 



where the 'magic' separable kernel K(t, t') is defined in [5] as 

2^nJ n {t)J n {t') + 2 

k=l 1=0 



K(t,t') 



+ C 2 k+l,2l+2{g)J2k(t)J2l+l(t') 



71=1 



(3.3) 



and its modification K*(t,t') as (13. 3p with the only replacement J\(t') — > Ji(t') — t'/2. In writing (13.21) 
the various contributions to the forcing term are separated according to their power of In s. The term 
proportional to Ins enters the BES equation for the cusp anomalous dimension, the four subsequent 
terms - independent of s - appear in the equation for the density which determines the virtual scaling 
function [21]. On the other hand, the contributions proportional to (lns) _n , n > 1, come from the term 
containing 



L-2 



P(s, g,t)=} j [cos tu h - 1] . 
h=l 

Expanding (13.41) for high s we have, because of (12. 5p . 

P(s,g,t) = Y i V n (g,t)(]ns)- n 



(3.4) 



(3.5) 



n=l 



Using again the Faa di Bruno formula, we can give a reasonably explicit expression for V n (g,t). After 
some calculations we obtain 

L-2 ^ i 

n ^ 2 lTm=l iP i m,h)' >rn n—r+1 n—r+1 

V n (g,t) = ^2t r COS y ^ rrn-r+1 ■ , ' Yl i™ = r ' Y mJm = U - ( 3 - 6 ) 

{j u ...,j n - r+1 } llm=l I™' m=l m=l 



r=l 



n n " rtl 7 ! 
llm=l Jm- 
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Remark 1 From the structure of the various a n ^ (see e.g. ( 12. 9p for n = 1, 2, 3) it follows that P(s, g, t) 



L-2 



is zero if S(m, L) = ^ (2/i + 1 — £) m = 0, Vm G N. This happens for L = 2,3. Thus, we conclude that 



h=l 



for L = 2, 3 the high spin expansion of the anomalous dimension does not contain inverse logarithmic 
powers in the spin, i.e. ^ n \g, 2) = ^ n \g, 3) = 0, Vn > 1. 

Remark 2 If we extend formally the result for S(m, L) to the case L = 1, we find again S(m, L = 1) = 0. 
We will use this formal property later, when discussing the relation to the large twist limit. 

It is now natural to introduce the high s expansion for the function S(k), 

oo 

S(k)= J2s {n \k)0ns)-\ (3.7) 

n=-l 

which allows to obtain a linear integral equation at each order in (lns) _n . Let us focus on n > 1. We 
immediately remark that S^(k) = 0, because the 1/lns term is absent in the large s expansion of 
P(s,g,t). The next step is the Neumann expansion for S^ n '(k), which is a standard procedure [7] in 
the case of an integral equation with separable kernel, 

S {n \k) =f^ S W(g)l&^ ^\g,L) = V2gS[%) . (3.8) 
P =i 



The Neumann expansion transforms the linear integral equation for S^ n \k) into a linear infinite syste 1 ^ 

- +0o «P„(g,i)^ p -i(v / 29() 



siUs) = -Wp-D / 7 s^Ta — -1>- 1 )EW»tf(t). 

2 m=l 



in) ^ ~ -2 P r ^MM^fft) -A P ±Z 2p , m (9)(-irS^(9), (3.10) 

2 m=l 

for each of the Neumann modes Sp n \g). A look at (13.101) shows that Sp n \g) are expressed in terms of 
the "reduced coefficients" S^ n \g) - defined in (4.23, 4.24) of |23 - and the quantities (J2J| a^'\ with 
n' < n — 3: this makes it possible to build up a recursive calculation scheme, opening the way to push 
the computation up to the desired order in Ins, in a way similar to [Mj. And, indeed, from (13.101) we 
can get the expression for ^ n \g) in terms of S[ n \g) and a H: h as 

L-2 n-r+1 

m,h)^ m n— r+1 n— r+1 



= -2nJ2s[ r/2 \g) cosy T-r+i ' Z Jm = r ' m ^ = n & lV > 

r=1 Uu-jn-r+l} Yl jm} m=1 m=1 



m=l 



6 We use the notation: 

7 f + °° dt J n (V2gt)J m (V2gt) 

Z n>m {g)= — . (3.9) 
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where a m ,h in (13.111) depend on \ with n' < n — 3. Explicitly, we have 7 ( - 1 )(g, L) — and 



7 (2) (^) 
7 (3) G?,£) 
7 (4) (^^) 



V2g o/ *_ 1)v> (L - 3)(L - 2)(L - 1) S^fo) 
o 



3(7, 



-2V2g y f° (L - 3)(L - 2)(L - 1) S?\g) 
3W ) 3 



\72#2 



7T 



7T 



(1) 



3 (a, 



(Oh 2 



3 («£">)» 2(^") 2 



(L-3)(t-2)(t-l)- 



7T <X 



90 K~^ 5 



7T 



360 (tr, 



2 j^(L - 3)(L - 2){L - 1)(5 + 3L(L - 4))j S^fo) 
(L - 3)(L - 2)(L - 1)(5 + 3L(L - 4)) S{ 2 \g) 



(-i) 
i 



(-1), 
o 



(3.12) 
(3.13) 



(3.14) 



These expressions are valid Vg and will be disentangled in the weak and strong coupling limit. Weak 
coupling expansions are provided in Appendix A, the strong coupling leading term in Section 5. 



4 Relation to the large twist case 

The expressions for 7 < - n )(g, L) can be written in terms of the 'generalised' scaling functions f n (g), fn (g), 
describing - according to formula (11.31) - the anomalous dimension in the limit fll .4jl . 

Beside the widely investigated f n (g), the functions fn\g) were studied in [26]. Expressions for 
fn\g), t > 1, can be obtained following the lines of Appendix A of [26] and are explicitly given, 
when 1 < r < 4 and 0<n<4 — r, in Appendix B of this letter. In the limit L — > oo expansion 
(ED coincide!] with f ll .3D . after using the definition j — [L — 2)/ Ins. Therefore, from the comparison 
between (JT72J) and (11.31) . one has, at large L, 



7 ^(g,L) = (L — 2) n+1 f n+1 (g) + £(L - 2) n ~ r fit(g) , n > 1. 

r=0 



(4.1) 



Since ^ n \g,L) is a polynomial in L — 2, (14.11) can be extended to arbitrary L. Now, we use the 
fact that , ~f^(g,2) = j^ n \g,S) = 0, (see Remark 1): this implies, when n > 1, f$ (g) = and 

fn+i{g) + X] fn-r(d) = O5 respectively. In addition, the 'formal' property announced in Remark 2, 



r=0 



7^(g, 1) = 0, produces the other relation (— l) n+1 / n+1 (g) + XX - 1)™ r fn-Ad) = 0, n > 1. Imposing 



r=0 



such relations on the general expression (14.11) . one gets 7 < - 1 ^(g, L) = and, Vn > 2, 

n— 3 



7 (n) (<7,£) 



(L-2) 



n+l 



(L-2) 3 - ± ^ i 



f n+1 (g)+J2 



r=0 



(L - 2) n ~ r -(L-2) 



3+(-l)" 



(4.2) 



7 This is true because we include all the terms (Ins) ™ in both expansions (the next term would be of order 0(lns/s)). 



S 



which expresses 7^(g, L) in terms of the generalised scaling functional at fixed j. 

Let us now show that (14.21) agrees with formula? (13.121 - I3TT41) . We first notice that entering the 
expression of "f^ig, L), is zero, since the density at the order (Ins) -1 vanishes. Then, using notations 
of [26] (cf. formulae (4.13,15) of that paper), we remark that = (L — 2)a^~ 1 ' 1 ^ + <j( 0,0 \ Plugging such 
expression into (13.121 - l3~l~4l) and using results from [26] and Appendix B (formula ( 1B.2j) ). we obtain 

^\g,L) = [(L-2) 3 -(L-2)]/ 3 (<?), (4.4) 

^ 3 \g,L) = [(L-2) 4 -(L-2) 2 ]/4(^) + p-2) 3 -(L-2)]/3 (0) (^), (4.5) 

^(g,L) = [(L-2) 5 -(L-2)]f 5 (g) + [(L-2Y-(L~2) 2 ]jt\g)+ (4.6) 

+ [(L-2)*-(L-2)]fi%), 

which indeed agree with (14. 2\\ . 



5 Strong coupling limit 

We are now in the position to extract the leading strong-coupling behaviour of (13. 12113.14"!) in an analytic 
way. To this purpose we briefly recall some results obtained in [21]. At strong coupling, the first 
component of the solution of the "reduced system" can be written as 

V2gS[ n \g) = { --^- ( - ) m(g) + 0(e"*), n>\ (5.1) 



2tt 



and, for the BES related densities a\ q , we get 



< 1} = -t (|) 2 " m(g) + O (e"^) , q > , (5.2) 



where m(g) is the mass gap of the 0(6) Non-Linear Sigma Model embedded [27] in M = 4 SYM 

m ^) = Y^9 1/4 [l + 0(g- 1 )]e-^ + ... . (5.3) 

Let us consider the remaining densities of Bethe roots. As far as 7*- 1 - ) (g, L), . . . , j^(g, L) are concerned, 
we only need to know <J^\ since is zero. Using results of [261 121], we get 

g _> oo 4 0) = -4\ng + 0(g°) . (5.4) 



8 The 'inverse' relation, which expresses the generalised scaling functions in terms ofj( n \g, L), comes from the deriva- 
tives of the latter with respect to L: 

f%- m) {9) = —^l {n) {9,L)\ L= z, n>l, 0<m<n+l. (4.3) 
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Hence, taking into account the previous equations, the leading strong coupling limit of (I3.12H3.14I) turns 
out to be 



7 (2) 

7 (4) G?,£) 

6 Summary and conclusion 



(L-2)] + 0(e-3t) , (5.5) 
2) 3 -(L-2)]+0(V^) , (5.6) 
(L -2)]+0 In g e^f ) . (5.7) 



7T 



24m(g) 
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3m 2 (g) 



[(L - 2f - 
hxg[{L- 



\n 2 g 

m 3 (g) 



[(L - 2) J 



In this letter we have studied the sub-logarithmic terms 7^(<?, L), n > 1, appearing in the high spin 
expansion ( II. 2p for the minimal anomalous dimension of twist operators in the sl(2) sector of A/" = 4 
SYM. We have found the general expression (13.111) . after solving via a recursive procedure the linear 
systems (13.101) : in particular, we have proved that ^ l \g,L) = and that ^ n \g,2) = ^^(g^S) = 0, 
Vn. Then, we have found the connection (14.2jl between j( n '(g,L) and the generalised scaling functions 
appearing in the limit ( 11.41) . Finally, weak and strong coupling limits for j^ n \g, L), with n — 2, . . . 4, 
are provided in Appendix A and in (15.51 - IBTTl) . respectively. In the large coupling regime and at leading 
order we would like to remark the appearance in ^ n \g,L) of the mass gap for the 0(6) Non-Linear 
Sigma Model: it would be nice to discover if a motivation for this result could be found on the string 
theory side of the correspondence, somehow miming [28] . 

As for the future, the connection between ^ n \g, L) and the generalised scaling functions describing 
the limit (jl.4p is worth to be investigated, especially at strong coupling. Another possible development 
is the study of the entire si (2) spectrum of anomalous dimensions [29], in the spirit of comparison to 
the recent developments on spiky strings [30] . 
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A Weak coupling limit 

The weak-coupling expansions for y( n \g,L), n = 2,3,4 are given in the present appendix up to the or- 
der g 8 . Such expansions can be easily obtained from the linear systems (|3.10p using a program of sym- 
bolic manipulation, in this case Mathematica® . With little computational effort we were able to reach order 
g 22 : for obvious reasons we do not give such higher orders here, but they are available in the web page 
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http:/ / www-fp.usc.es/~grinza/gamma/gamma.html 



^ 2 \g,L) = [(L-2f-(L-2)] 
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B Computations of nonlinear terms 

In this appendix we report on the computations of the generalised scaling functions fn\g), with 1 < r < 4 

and < n < 4 — r. In order to do that, we have to start from equation (13. 2D . Forcing terms which contribute 

L-2 

to those generalised scaling functions are the ones containing [ cos t u h ~ 1] • We evaluate such an expression 

h=l 

in the limit (jl.4p by using results contained in Appendix A of [26]. Using formulas (A.3) and (A. 5) of [26], we 
first have that 
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where cr m (c) denotes the m-th derivative of the density cr(v) in v = c. Restricting to the cases of our interest: 
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where we used notations of [26]. After writing the systems for the Neumann modes we eventually realise that 
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